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ABSTRACT 

Aims. To study the usefulness of CF + as a tracer of the regions where C + and H 2 coexist in the interstellar medium. 

Methods. We used the Plateau de Bure Interferometer to synthesize CF + J= 1-0 absorption at 102.6 GHz toward the core of the distant 
HII region W49N at 1 = 43.2°, b=0.0°, and we modeled the fluorine chemistry in diffuse/translucent molecular gas. 

Results. We detected CF + absorption over a broad range of velocity showing that CF + is widespread in the H 2 -bearing Galactic disk 
gas. 

Conclusions. Originally detected in dense gas in the Orion Bar and Horsehead PDR, CF + was subsequently detected in absorption 
from diffuse and translucent clouds seen toward BL Lac and 3C111. Here we showed that CF + is distributed throughout the diffuse 
and translucent molecular disk gas with N(CF + )/N(H 2 ) = 1.5 - 2.0 x 1(T 10 , increasing to N(CF + )/N(H 2 ) = 3.5 x 1(T 10 in one cloud 
at 39 kms -1 having higher N(H 2 ) « 3 x 10 21 cirT 2 . Models of the fluorine chemistry reproduce the observed column densities and 
relative abundance of HF, from which CF + forms, but generally overpredict the the column density of CF + by factors of 1.4-4. We 
show that a free space photodissociation rate T ^ 10 -9 s -1 , comparable to that of CH, might account for much of the discrepancy but a 
recent calculation finds a value about ten times smaller. In the heavily blended and kinematically complex spectra seen toward W49, 


CF + absorption primarily traces the peaks of the H 2 distribution. 
Key words, interstellar medium - abundances 


1. Introduction 

Recent advances in mm- and sub-mm astronomical spectroscopy 
open the possibility of tracing Hi in interstellar gas over a 
much wider range of column density and Hj-fraction, even when 
carbon monoxide, the usual surrog ate tracer of H? is not ob¬ 
served. One new Ho-tracer is HF d Sonnent rucker et all 1201 ft 
iGerin et~aIll2012tlSonnentrucker et al.ll2015l) which is~predicted 
to contain much or m ost of the gas-phase fluori ne even when the 
Ho-fraction is small dNeufeld & WolfireLl2009l) . The abundance 
of flu orine in diffuse neutral gas is well-established (ISnow et al.L 
120071 ) and the HF/H 2 ratio is relatively insensitive to physical 
conditions, so that N(HF) is a useful estimator of N(H 2 ). 

Ground-state absorption from HF is observed at 1.2 THz. 
Observing this is a challenge now, but HF forms the ba¬ 
sis for an observable fluo rine chemistry that i s accessible at 
mm-w avelengths via CF + dNeufeld et all 120061 : iGuzman et all 
l2012allbh whose ground state J=l-0 rotation transition is at 102.6 
GHz. CF + is an immediate descendant of HF, formed in the re¬ 
action C + + H F —» CF + + H and destroy ed by recombination 
with electrons dNeufeld & Wo I fire) l2009t) . Because C + is the 
dominant form of gas-phase carbon well beyond Ay = 1 mag 
dSofia et al.L 120041 : iBurgh et al.Ll2010t) , and because the CF + /HF 
and CF'/fL ratios are easily modelled, CF + may also be a use- 
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ful tracer of molecular gas where C + is abundant. It thus has the 
potential to identify the fraction of C + arising in molecular gas. 

C F + was discovered in a targeted observation of the Orion 
Bar dNeufeld et all 120061) , subsequentl y detected in a 43mm 
spectral sweep of the Horsehead PDR dGuzman et al.L l2012ah 
and then detected in a 43mm spectral sweep of the absorption 
spectra of diffuse and tran s lucent clouds occult i ng BL Lac and 
3C111 dLiszt et al.L 120141) . iFechtenbaum et all (12015l) recently 
detected CF + emission toward a high-mass protostar in a similar 
spectral sweep and used it to infer a source of carbon ionization 
in the collapsing protostellar envelope. 

Here we show that CF + is observable in gas seen along the 
line of sight to the giant HII region W49 N at / = 43.16°, b = 
0.02° at a distance of 11.1 kpc ( Zhang et all 2013h and more 
generally in the gas that gathers near the terminal velocity at v ^ 
60 kms -1 due to velocity-crowding ( Burton. 1971 ). As a result, 
CF + is clearly established as a commonly-distributed species in 
diffuse molecular gas. 

The structure of this work is as follows. In Sect. 2 we de¬ 
scribe the new and existing observational material that is dis¬ 
cussed. In Sect. 3 we present the detection spectrum of CF + to¬ 
ward W49N and compare its kinematics and column densities 
with those of other species seen in absorption along the same 
sightline. In Sect. 4 we compare the column densities of HF and 
CF + with models of the fluorine chemistry in diffuse/translucent 
clouds and Sect. 5 is a summary. 
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Fig. 1 . Line profiles observed toward W49N. Top: emission 
spectra from the ARO 12m antenna at « 1' spatial resolu¬ 
tion. Brighter emission below 20 km s -1 from the W49N H II 
region-molecular cloud complex has been truncated to allow bet¬ 
ter viewing of weaker emission from the Galactic plane clouds. 
Middle: absorption from weaker-lined transitions including CF + 
studied at the PdBI. Bottom: absorption in stronger and more 
ubiquitously distributed species observed at the 30m. The ab¬ 
sorption data are plotted as line/continuum ratios, displaced as 
needed to fit the plot. In some cases the profiles have been scaled 
to make weak features more easily visible, with a scale factor as 
shown at the right hand side of the spectrum. 


2. Observations and data reduction 

2. 1. PdBI absorption measurements toward W49N 

We observed the CF + 1=1-0 line at 102.587 GHz (Table 1) with 
the Plateau de Bure Interferometer during 34 hours of telescope 
time, corresponding to 9.3 hours of on-source time scaled to 
a six antenna array, after filtering out low-quality visibilities. 
The observations were carried out in August and September 
2013 in the D configuration, achieving a spatial resolution of 
5.95"x 4.72" so that 1 Jy/beam is equivalent to 4.14 K. The 


typical precipitable water vapor amounted to 6 mm with a typ¬ 
ical system temperature of 150 K. The bright quasars 3C454.3 
and 3C279 were used to calibrate the bandpass. The phase and 
amplitude temporal variations were calibrated by regularly ob¬ 
serving two nearby quasars (1827+062 and 1923+210). The ab¬ 
solute flux scale was derived from observations of MWC349. 
The calibration and imaging of the PdBI data was done with the 
GILDAS0/CLIC and MAPPING software. The data were gridded 
onto a 1.24"pixel spatial map with rms noise 0.0105 Jy/beam in 
spectral channels of width, spacing and resolution 78.1 kHz or 
0.228 kms~'. 

The 3mm continuum peaked at a(J2000) = 19*10™ 13.45', 
d>(J2000) = 9°06'13.2" or 1 = 43.2°, b = 0.0° in Galactic co¬ 
ordinates, with a flux of 3.59 Jy/beam. The CF + absorption line 
profile was extracted at a position with a cleaner spectral base¬ 
line that was displaced one pixel to the West and South hav¬ 
ing a flux of 3.47 Jy/beam, leading to a channel-channel rms 
line/continuum ratio of 2.4 10 3 . The CF + spectrum is shown 
in the middle panel of Fig. 1. The CF + absorption is overlain 
by the emission profile of an H66e recombination line from the 
ionized gas and the procedure for extracting the CF + absorption 
profile in the presence of the recombination line is described in 
Appendix A. 


2.2. Other absorption profiles toward W49N 


For comparison with other species observed in the 3mm band 
we used the absorption profilcs_of HNC and HCO + taken at 
the IRAM 30m telescope by Godardjhal.j 201 Til) , along with 
similar profiles of ^h2, CtH ( Gerinet all 2011 ) and H l3 CO + . 
These spectra have 39 kHz channel spacing and 78 kHz spec¬ 
tral resolution equivalent to 0.117 kms -1 spacin g at 100 GHz. 
The HCO spectrum shown in Fig. 1 is from lLiszt et al.i (20141) . 
shown at 0.270 kms -1 channel spacing and resolution. For CH 
we used the hyperfine-deconvolved 532 GHz absorption profile 
from iGerin et al.l (1201 Ol) at 0J_35 km s_|_ spacing and resolution, 
for HF the profile of Sonnentrucker et aTI (120101) at 0.141 kms~' 
spacing and resolution and for C + a smoothed absorption spec¬ 
trum from iGerin et al.l (2015 ) at 0.63 km s 1 spacing and resolu¬ 
tion. The sub-mm profiles were taken as as part of the Herschel 
PRISMAS Project. The /13mm profiles are shown in Fig. 1 and 
with the spectrum of CH in Fig. 2. Profiles of the C + and HF ab¬ 
sorption lines can be seen in Fig. 4. Line frequencies and optical 
depth-column density conversions for these species are given in 
Table 3. 

When profiles are compared on a channel-by-channel basis 
the narrower of the two was convolved to the spectral resolution 
of the other using a gaussian kernel and then regridded using 
cubic splines to create a common velocity axis. 


2.3. CO and HCO + J=1-0 emission 

The emission profiles shown in the top panel of Fig. 1 were 
acquired at the original ARO 12m telescope during April of 
2011 by position switching against a distant off position whose 
spectrum was determined by position switching against another 
similarly distant off-position and added back to the data. The 
beamwidth of this antenna was 1' at 110 GHz, scaling inversely 
with frequency. The temperature scale of the ARO 12m tele¬ 
scope, like that of the NR AO 12m telescope before it, was T*, 
approximately 0.85 T m b. The emission spectra have spectral res- 

1 See http://www.iram.fr/IRAMFR/GILDAS for more informa¬ 
tion about the GILDAS software i l*e1\ . 2005 ). 
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Fig. 2. Histograms of integrated optical depth over discrete ve¬ 
locity ranges for species observed in absorption toward W49N, 
normalized to a maximum of unity in each case. The quantities 
plotted are those given in Table 2. Scaled and offset absorption 
spectra are shown superposed in red to illustrate the underlying 
kinematic features in the velocity bins. 


olution and channel spacing 48.8 kHz. The HCO + spectrum 
shown in Fig. 1 is the average of four spectra taken at posi¬ 
tions displaced +1.4' to the East and North, to avoid contami¬ 
nation of the very weak emission by absorption. The CO pro¬ 
files were taken directly toward the continuum, as the W49N 
continuum was a small fraction of a Kelvin at the ARO 12m 
antenna. Reliably observing carbon monoxide in absorption to¬ 
ward W49N would require spatial resolution and resolved flux 
substantially exceeding those achieved in the PdBI synthesis of 
CF + here, for which the peak flux of 3.6 Jy/beam was equivalent 
to 15 K. This is not so much larger than the 4 K emission bright¬ 
ness of CO that a pure absorption spectrum could be acquired. 


2.4. Optical depths, error estimates, upper limits and column 
density conversions 

The spectroscopic parameters of the observed CF + transitions 
are given in Table 1 and optical depth integrals for CF + and other 
species discussed here are given in Table 2. Quantities in paren¬ 
theses in Table 2 are l<x rms errors in the tabulated quantities 
determined from the rms noise level at zero absorption but ac¬ 
counting for the saturation in the absorption line profile. Where 
the integrated optical depth was detected at a level well below 
3 <x the integral is shown as zero and 3<x upper limits are plotted 
in the various Figures. 

In the present experiment the hyperfine structure of CF + was 
unresolved and the integrated optical depth is 50% larger than if 
the stronger component had been observed alone. Hence N(CF + ) 
is derived from the integrated optical depth in Table 2 using 
N(CF + ) = 2/3 x 2.02 x 10 13 cm' 2 f rdv with the optical depth 
integral expressed in km s' 1 . Conversion factors between inte¬ 
grated optical depth and column density are given separately in 
Table 3, assuming excitation in equilibrium with the cosmic mi¬ 
crowave background. 


3. The line of sight toward W49N 

3. 1. Line profiles 


/13mm line profiles observed toward W49N are shown in Fig. 1. 
Emission lines from dense gas in the W49N complex peak at 0 
- 15 kms' 1 but may have broad wings extending well beyond 
this range: CF + and HCO are notable exceptions because they 
do not appear in emission at the PdBI. When the emission is not 
too broad the continuum and the source emission profiles are 
noticeably absorbed at 5-20 km s' 1 . Much of this absorption is 
due to gas within and on the near-side of the W49N complex 
itself: The relationship between emission and absorption pro¬ 
files of HCO and H 13 C O + as se en by the PdBI toward W49N 
is shown in Fig. 4 of iLiszt et al.l (120141) . Interstellar absorption 
from gas lying closer to the Sun must also be present at v $ 20 
km s' 1 but is masked by the contribution from W49N. As shown 
in Appendix B, the velocity range v < 20 km s' 1 corresponds to 
the nearest 1.75 kpc along the line of sight. 

At bottom in Fig. 1 are strong interstellar absorption lines 
of HNC, HCO + and 3h2, representing three distinct chemical 
fami lies. HCO + ( Lucas & Lisztl 119961) and the hydrocarbons 
CH dGerin etakl boiQl) . C?H and ?h2 (iLucas & Lisztl 120001: 
iGerin et al.Ll2011l) are especially ubiquitous and show the widest 
range of absorption. The hydrocarbons and HCO + have rather 
stable abundances with respect to H 2 , which we use to derive col¬ 
umn densities of Ho using X(CH) = N(CH)/N(H 2 ) = 3.5 x 10' 8 
(tSheffer et all i2008t) to ward W49N and X(HCO + ) = 3 x 10' 9 
(Table 4 and ILiszt et al.l (l2010t) ) for quasar lines of sight where 
CH was not observed. Scatter in the CH-H 2 relationship, uncor¬ 
rected for measurement errors, is a factor 1.6 as employed in Fig. 
5. The column densities of CN, HNC and HCN vary together in 
fixed proportion but only achieve high abundance relative t o H 2 
at higher N(H 2 ) (iLiszt & Lucasl I200U Godard et al l 201 Oil and 
presumably at somewhat higher number density as well. For this 
reason they have more limited absorption profiles but serve to 
distinguish the more diffuse molecular gas. 

In the middle panel of Fig. 1 are profiles of three weaker- 
lined species including the detection spectrum of CF + . With 
sufficient sensitivity the interstellar H 13 CO + and HCO + absorp¬ 
tion profiles would presumably be identical in shape. With the 
present sensitivity all of these species show absorption over only 
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Fig. 3. Column density of CF + vs. those of HCO + (left), C 2 H (middle) and HNC (right) for the velocity intervals in Tables 2 and 4. 
Shown for comparison are data from diffuse cloud sightlines toward BL Lac and 3C111, and for the Horsehead (TdC). The small 
brown points are the individual channel values of optical depth in the original spectra, scaled equally in each coordinate to fit in the 
plot area. Upper limits are 3 cr. The solid lines are regression lines fit to the detections toward W49, ignoring the upper limits. 


a limited range, most notably in the 40 km s 1 cloud feature, but 
also in the terminal-velocity ridge above 60 kms" 1 . 

At top in Fig. 1 are J=l-0 emission profiles of CO ( 12 CO), 
13 CO and HCO + . It is notable that CO emission is detected at all 
velocities at which absorption is observed in the hydrocarbons 
and HCO + . 13 CO is detected over a more limited range, in part 
due to considerations of signal/noise for the weaker line. The 
emission at 33 km s' 1 has a very high ratio of CO to 13 CO emis¬ 
sion, presumably reflecting optically thin emission and a small 
CO column density. 

As shown in the uppermost panel of Fig. 1, HCO + emission 
is detected at 39 km s 1 and v ^ 50 km s -1 at a level ~ 1 -2% that 
of CO as is typic al of emission observed in the gala ctic plane in 
the inner Galaxy dLiszti 119951 : iHelfer & Blitd . ll997f) . 


3.2. Line of sight kinematics 

The kinematic characteristics of the absorption profiles toward 
W49N result from several influences as illustrated in Fig. B.l 
and discussed in Appendix B. W49N is 11.1+0.8 kpc distant 
from the Sun (IZhang et ah , 1201 3l) . hence at galactocentric ra¬ 
dius R = 7.6 kpc, just inside the Solar Circle on the other side 
of the Galaxy for a Sun-center distance of 8.5 kpc. Hence lines 
originating within the W49N complex are closest to 0-velocity 
because of its inadvertent proximity to the Solar Circle on the 
far side of the Galaxy. The highest velocities appearing in emis¬ 
sion or absorption in Fig. 1 correspond quite well to the maxi¬ 
mum velocity predicted for a flat rotation curve with Rq = 8.5, 
&(R) = 220 kms -1 as shown in Fig. B.l. An apparent concen¬ 
tration of gas arises at the high-vel ocity end of the line profiles 
as the result of velocity-crowding ( iBurtonJ 197 ll) . as the slope of 
the velocity-distance relationship approaches 0 at the so-called 
sub-central point where the line of sight passes nearest to the 
center of the Galaxy. Gas is concentrated into the velocity range 
30 - 40 kms -1 , with notable gaps at higher and lower velocity, 
owing to sharp changes in the velocity-distance relationship re¬ 
sulting from spiral-arm streaming. This has not been extensively 
modeled for W49N but a detailed discussion of the influence of 


spiral arm streaming along the sightline toward W43 is given in 
iLiszt et all(ll993l) . 

3.3. Column densities and cloud properties 

To give a broad characterization of the properties of the gas ly¬ 
ing along the sightline to W49N we divided the velocity range 
of the interstellar absorption into a set of intervals based on the 
apparent features in the absorption profiles, and we calculated 
the integrated absorption optical depth for a variety of molecules 
shown in Fig. 1 and 2, as given in Table 2. It is tempting to iden¬ 
tify narrow-lined spectral features with individual clouds, espe¬ 
cially at 39 kms -1 , but there is severe blending at all velocities. 
The species in Table 2 include the three hydrocarbonsCH, CLH 
and 3h2 that are known to track each other dGerin et al.l 1201 ll) 
and to trace H 2 , and HCO + that also has a relatively stable rela¬ 
tive abundance with respect to H 2 . The hydrocarbons and HCO + 
have much more broadly distributed absorption than HNC, for 
instance, even though the peak absorption is not much different. 

The results are shown graphically in Fig. 2 where the his¬ 
tograms of integrated optical depth and the line profiles are 
shown superposed. In each case the histogram is normalized to 
have a unit peak. Species which most closely track the total col¬ 
umn density of H 2 without regard to the number or column den¬ 
sities of individual regions (the hydrocarbons and HCO + ) have 
higher column density over the terminal velocity ridge around 
v = 60 kms . The hydrocarbons and HCO + also have a rela¬ 
tively high ratio of absorption at 33 km s 1 compared to that at 39 
km s 1 . HNC, which is known to be more abundant in individual 
clouds of higher H 2 column density dLiszt & LucasL l200l[) . be¬ 
haves oppositely to the hydrocarbons and HCO + . HCO is an ex¬ 
treme example of a species appearing preferentially in the more 
diffuse gas in the terminal velocity ridge. CF + is clearly most 
similar to HNC, which is something of a surprise given its origin 
in HF, which tracks H 2 , and C + that is ubiquitous (see Sect. 4 
and Fig. 4). 

The optical depth integrals in Table 2 are converted to rela¬ 
tive abundance N(X)/N(H 2 ) in Table 4 using the column density- 
optical depth conversion factors given in Table 3 to form the ra- 
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tio N(X)/N(CH) and N(CH)/N(H 2 ) = 3.5 x lO^ 8 . Shown at the 
bottom of Table 4 are two summary quantities: the mean rel¬ 
ative abundance found by integrating over the profiles of the 
species listed relative to the total CH column density, and the 
mean and rms of the individual entries in the Table, for those 
int ervals where the species in question was detected. As noted 
bv iGerin et al.l (12011 1 the hydrocarbons have quite stable abun¬ 
dances relative to each other: CL H and 3h2 have rms of 20% 
or less of their mean table entries. HF is also relatively stable 
with respect to CH, with a fractional rms of 23%. The fractional 
rms of HCO + and HCO are intermediate, 42%, and those of CF + 
and HNC are the largest, 48% and 69%, respectively. The over¬ 
all average abundance X(HCO + ) = 2.7 10~ 9 is consistent with 
the value X(HCCV_) = 3 10 9 derived from quite other consider¬ 
ations bv [Liszt et al.l (l2010h . 

In all, the results in Table 4 shows an impressive consistency 
of the diffuse cloud chemistry over the Galactic disk, and val¬ 
idates the empirical methods that have been used across many 
spectral domains to extract molecular column densities and rel¬ 
ative abundances. 



Fig. 4. Optical depth profiles of CF + (scaled upward by a factor 
50) and the species HF and C + that form it. 


4. The abundance and chemistry of CF + 

4.1. N(CF + ) and chemical relationships 

Extending the analysis in [Liszt et aTI (120141) . Fig. 3 shows the 
CF + column density with respect to HCO + , C 2 H and HNC for 
W49N and for the sightlines observed toward compact extra- 
galactic continuum sources BL Lac and 3C111, as well as data 
taken toward the Horsehead PDR. Shown for W49N are the in¬ 
tegrals tabulated in Table 2 and the cloud of channel-by-channel 
datapoints equally scaled in both dimensions so as to overlay the 
datapoints while preserving the slope. Shown in each panel is a 
regression line fit to the W49N detections alone. The slopes of 
these fits are 0.78 ± 0.56 for HCO + , 1.97 ± 0.44 for C 2 H and 
0.70 + 21 for HNC. 

The velocity ranges chosen toward W49N have column den¬ 
sities of HCO + , C 2 H and HNC overlapping those seen toward 
BL Lac and 3C111 and the regression line fits show a high de¬ 
gree of consistency between the W49N sightline, 3C111 and, 
where possible, the Horsehead PDR. The line of sight toward 
BL Lac is an obvious outlier, with 2-3 times more CF + compared 
to the other datapoints. The 39 kms -1 cloud toward W49N is a 
lesser outlier with respect to HCO + . 

A striking aspect of Fig 3 is the quadratic slope of the CF + - 
C 2 H relationship as compared with the much slower than linear 
variation with respect to HNC. This is indicative of a rapid vari¬ 
ation o f N( HNC) with N(C 2 H), as was in fact shown in Fig. 3 
of lLiszt & Lucas! ([200 ll) : N(HNC) increases abruptly by about a 
factor 100 in individual cloud features for N(C 2 H) L 10° cm~ 2 , 
and also for N(HCO + ) L 10 12 cm -2 . By contrast, CF + has a faster 
than linear variation only with C 2 H. 

4.2. Modelling the fluorine chemistry 

The fluorine ch emistry of diffuse a n d tran slucent gas has been 
summarized by iNeufeld & Wo I fire! ( 2009 ). It is dominated by 
the formation of HF, which is expected to contain much or most 
of the free gas-phase fluorine in reg ions where the H 2 fraction i s 
appreciable. However, as noted by ISonnentrucker et ahl (12015l) . 
the rate constant for the primary reaction governing the forma¬ 
tion of HF, F + H 2 —» HF + H, has recently been found to 
be rather sm aller at low temperature than previously believed 
dTizniti et alll2014 ). This complicates discussion of the fluorine 


Table 1. CF + transitions observed and column density-optical 
depth conversions'* 


Frequency 

MHz 

F' - F 

J'-J 

A fi,f 

10 6 s -1 

N(X)/J Tdv b 
cm~ 2 (km s-V 

102587.189 

1/2-1/2 

1-0 

4.82 

4.04 x 10 13 

102587.533 

3/2-1/2 

1-0 

4.82 

2.02 x 10 13 


° Spectroscopic data from [Guzman et all i20T2bj) and 

www.splatalogue.net 

b Assuming excitation in equilibrium with the cosmic microwave 
background 


chemistry because the HF/H 2 rat io is more variable, and ge ner¬ 
ally smaller, than in the mo dels of lNeufeld & Wolfir3 (120091) . see 
ISonnentrucker et akl (l2015h . 

HF is also formed in a minor way by the reaction of neu¬ 
tral fluorine atoms with CH, OH, and, more slowly, with wa¬ 
ter. In diffuse gas, HF is destroyed in reaction with C + (the 
primary process, forming CF + ), in reaction with He + (created 
by cosmic-ray ionization) and by photodissociation. CF + is de- 
st royed only by elect r on rec ombination in the chemical model 
of lNeufeld & Wolfirel d2009t) that was adopted here. Because C + 
usually dominates both the formation and destruction of CF + 
in diffuse gas, the CF + /HF ratio is generally insensitive to the 
number density as long as C + is the dominant form of gas-phase 
carbon and provides the majority of free electrons. 

The species involved in the CF + formation chemistry are ob¬ 
servable and Fig. 4 shows the optical depth profiles (equivalently 
the profiles of column density or column density per unit veloc¬ 
ity) of C + , HF and CF + . Because it is so widespread and so heav¬ 
ily blended, and to some extent because it was observed with 
some three times lower velocity resolution, the C + profile varies 
little across the peaks in HF and CF + . The optical depth of CF + 
varies more strongly than that HF or the product r(C + )*r(HF). 

To model the CF + chemistry we adopted a_very slightly sim- 
plified version of the fluorine chemistry of Neufeld & Wolfire 
( i2009t) using the new HF formation rate of [Tizniti et al. ( 20141) 
and fixed gas-phas e elemental abundances [ Cl/fHl — 1.4 x 10 4 
( Sofia et al l EZQOj) . [F]/[H] = 1.6 X 10~ 8 ( Snow et al .. l2007h . 
We held fixed the empirically determined molecular abundances 
X(CH) = 3.5 x 10- 8 , X(OH) = 10~ 7 (IWeselak et all 1201 Oh and 
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Fig. 5. Model and observed CF + and HF abundances. Upper: the 
HF/FF ratio in models of HF and Hi formation at total densi¬ 
ties n(H) = 64, 128 and 256 cm -3 . Lower: observed and model 
HF and CF + column densities. The W49N data were placed in 
the plane using N(CH)/N(H 2 ) = 3.5 x 10~ 8 and the data for 
BL Lac and 3C 111 were positioned relative to HCO + assum¬ 
ing N(HCO + )/N(Hi) = 3x 1 CL 9 . Profile sums over the velocity 
intervals in Table 2 are shown as solid symbols and upper limits 
are 3 cr. Individual channel values of the HF and CF + optical 
depth spectra are shown as red and blue points, respectively, af¬ 
ter scaling in each coordinate by the Galactic line of sight veloc¬ 
ity gradient for a flat rotation cu rve (see text). The dashed light 
brown curve shows the results of lNeufeld & Wolfirel d2Q09t) . 


the H-Hi transition across 128 equi-spaced radial shells. X(HF) 
» 1.1 x 1CL 8 , independent of n(H) for N(H 2 ) S 2x 10 21 or Ay 
F 2 mag. At higher N(H 2 ) carbon recombines to neutral carbon 
and CO and the photodissociation rate of HF decreases, allowing 
X(HF) to increase to the point where HF cons ume s all of the free 
gas-ph ase fluorine (see also the models of ISonnentrucker et al.l 
(20151)). X(CF + ) remains nearly constant as the decrease in C + 
is compensated by falling electron density. The HF/CF + ratio 
increases, which we have traced in the models to falling temper¬ 
ature and faster CF + recombination. 

Shown for comparison is the earlier prediction for N(CF + ) 
of lNeufeld & Wolfira d2009l) . which includes progressive freeze- 
out of the gas-phase fluorine. Agreement of the two predic¬ 
tions for CF + is surprisingly good considering the differences 
in the underlying HF che mistry. The HF column densities in 
iNeufeld & Wolfira (i2009i) are about 50% higher than here at 
N(H 2 ) < 2 x 10 21 cm- 2 . 


4.3. Model and observed results for HF 


Shown in Fig. 5 are HF and CF + column densities over the 
velocity intervals tabulated in Table 2 and 4, and the clouds 
of channel by channel optical depths scaled to fit in the plot 
area. As noted in Sect. 2, N(H 2 ) is approximated as NHL) = 
N(CH)/3.5 x lO -8 , here shown with an assigned error in N(H 2 ) 
of a factor 1.6 corresponding to the scatter in the empirical CH- 
Hi correlation of lSheffer et al. ( 2008) . Also shown in Fig. 5 are 
our measurement s of CF + , along with pri or measurements in the 
Horsehead PDR (Guzman et all 20 12ah and in absorption to¬ 
ward BL Lac and 3C111 (Liszt et all 12014) . The quasar absorp¬ 
tion datapoints were placed horizontally by a ssumi ng a relative 
abundance N(HCO + )/N(H 2 ) = 3 x 10- 9 (see lLiszt et aH (12010) 
and Table 4). 

The results for HF are well explained by the models, and 
have a nearly constant abundance relative to Hi at N(H 2 ) F 
2 x 10 2 1 cnT 2 , X(HF) = 1.1 x 10 / 8 in agreement with the re - 


sults of llndriolo et al. (12013) and ISonnentrucker et al.l (12015) . 
As sampled in HF, whose abundance with respect to H 2 remains 
steady up to the highest inferred values of N(H 2 ), freeze-out of 
the fluorine and/or collapse of the carbon ionization to CO and 
neutral carbon must occur at values of N(H 2 ) beyond those en¬ 
countered along the line of sight to W49N. The absorption line 
of HF becomes heav i ly satu rated at higher column densities, see 
ISonnentrucker et aD (12010) 


4.4. Model and observed results for CF + 


X(HiO) = 1.2 x 10 8 (ISonnentrucker et all 12015) and embed¬ 
ded this chemistry in models of the heating and cooling of uni¬ 
form density gas spheres, such as were recently used to calcu¬ 
late abundances of Hi, HD and CO, etc. (Oszi. I2007L l20l5h : 
these use standard values of the radiation field and photodisso- 
cation rates, et c. We adopted a cons t ant ga s-phase abundance of 
fluorine as in (ISonnentrucker et all 12015) and did not include 
a progressive freeze-out of gas-phase fluorine, which apparently 
occurs only at higher values of the hydrogen column density than 
are encountered in this work, see Fig. 5 and the discussion below. 

Model results for HF and CF + are shown in Fig. 5 for three 
densities n(H) = 64, 128 and 256 cm -3 : the proportions of 
atomic and molecular hydrogen are calculated inside the mod¬ 
els self-consistently. The results plotted are integrated over the 
central line of sight toward models with varying N(H) following 


The predicted value over the range of N(H 2 ) observed toward 
W49 is X(CF + ) = 5.3 X 10 10 « X(HF)/20. This can be under¬ 
stood simply by comparing the rate of CF + formation from C + + 
HF —> CF + + H and the rate of CF + recombinat ion in a gas at 75 
K (the mean temperature of H 2 (kachford et alil2002t) ') in which 
C + provides 7/8 of the ambient electrons as determined from de¬ 
tailed examination of the internal structure of the models. 

The ratios N(CF + )/N(H 2 ) observed in absorption toward 
W49N (Table 4) are 3.3 x 10~ 10 for the gas around 39 kms’ 1 
and 1.40 + 0.16 x 10 10 for the other four velocity intervals in 
which CF + was detected. The smaller of the observed relative 
abundances is about a factor four below the model prediction 
while that at 39 km s _1 is about 40% too small. The line of sight 
toward BL Lac was already cited as an outlier in the discussion 
of Fig. 3. The CF + abundance along the line of sight through the 
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translucent gas toward 3011 is consistent with the smaller of 
the CF + abundances toward W49N. 

Agreement of the model and observed results for HF implies 
that the ratio of HF formation and destruction rates is correct 
under the assumption of equilibrium of the chemical rate equa¬ 
tions. Given that the f ormatio n rate using the recently measured 
reaction rates of lTizniti et alJ ( 2014 1 should be secure and the H- 
H 2 transition is modelled in detail, it follows that the destruction 
rate of HF should also be correct. Therefore the formation rate 
of CF + is also correct to the extent that formation of CF + dom¬ 
inates the destruction of HF. This implies that small CF + abun¬ 
dances arise from the absence of some important CF + destruc¬ 
tion mechanism in the models, and perhaps that this mechanism 
is suppressed in the 39 km s 1 gas toward W49N and toward BL 
Lac. 


4.5. A role for CF + photodissociation? 


We considered that this could point to a role for photodisso¬ 
ciation, which was not incorporated in the underlying chem- 
ical_model given the rapidity with which CF + recombines. 
[Guzman et al.l (12012al) noted that even a relatively large free- 
space photodissociation rate F = 10~ 9 s _1 was too slow to affect 
the abundances seen in the Horsehead nebula and they neglected 
it in their calculation of the CF + column density. 

To clarify this matter we examined our mod¬ 
els to see what were the actual rates of photodisso¬ 
ciation of CF + . The ranges of recombination rates 

in the models averaged along the line of sight were 


1.8 10 ..0.8 10“ 9 s _1 (64 cm -3 ), 3.6 10 9 ..1.5 lO^V 1 (128 cm” 3 ) 
and 6.9 10~ 9 ..1.9 10 _9 s _1 (256 cm -3 ) with the higher values at 
the smallest N(H 2 ) and vice versa. Thus a free space photodis¬ 
sociation rate T ^ 10 9 s“' would have a very noticeable effect 
on X(CF + ), introducing number and column density depen¬ 
dences beyond the very small ones in the models shown here. 
This is a r elatively high rate, but not higher than that for CH 
( iVan Dishoeck & Black . 19861) which is quoted in the UFDA 
database dMcElrov et al. ~ Eoi3l) as 1.7 x 10~ 9 s _1 However, 
we subsequently learned that the free-space photodissociation 
rate of CF + has recently been calculated to be 2 x 10 l() s 1 
(Dayou and Roueff, in preparation), effectively eliminating 
photodissociation as a practical solution to the discrepancy 
between the calculated and observed CF + /HF ratio. 

In closing we note that some uncertainties remain regarding 
the underlying chemistry. The adopted CF + recombination rate 
was measured in a gas of rotationally-warm CF + , and may need 
refinement. Moreover, the assumed formation rate for CF+ by 
reaction of C+ with HF has not yet been measured in the lab¬ 
oratory; here, we simply adopted the capture rate computed by 
NW09, assuming - in effect - that every capture results in a reac¬ 
tion. In reality, the CF+ formation could be significantly smaller 
than what we have assumed, providing an alternative explanation 
for the lower-than-predicted CF+ abundances reported here. 


5. Summary 

CF + , formed in the reaction C + + HF —> CF + + H, was orig¬ 
inally detected in the Orion Bar and subsequently detected in 
the Horsehead PDR and then in diffuse/translucent molecular 
clouds occulting BL Lac and 3C111 and in a massive proto- 
stellar envelope. In this work we detected the 102.6 GHz .1=1- 
0 lines of CF + (Table 1) in absorption from gas lying in the 
galactic plane toward W49N , a giant HII region at a distance 


of 11.1 kpc from the Sun at l - 43.2°, b = 0.0°. As detailed 
in Appendix A the CF + profile was extracted in the presence 
of an overlying H66e recombination line. As shown in Fig. 1 
and 2 and summarized in Tables 2 and 4, narrow absorption 
was seen in features at 33 and 40 kins' 1 corresponding to well- 
known discrete diffuse/translucent molecular gas clouds lying in 
the Galactic plane. But absorption from more diffuse material 
was also observed at intermediate velocities and in the terminal- 
velocity ridge at v ^ 60 kms -1 caused by velocity crowding 
of unrelated gas parcels as the line of sight velocity gradient 
approaches zero at the sub-central point (see Appendix B and 
Fig. B.l). This establishes CF + as a fairly commonly-distributed 
species in diffuse and translucent molecular gas. 

In Fig. 2 (see Table 2) we broadly compared the velocity dis¬ 
tribution of CF + with those of several hydrocarbons and HCO + 
tracing the total column density, with HNC that appears strongly 
only when the number and column densities are higher locally, 
and with HCO, an interface species that is known not to pen¬ 
etrate dense gas. The hydrocarbons and HCO + have higher in¬ 
tegrated opacity on the high-velocity side of the profiles where 
column density aggregates as the result of velocity-crowding in 
portions of the line of sight at the smallest galactocentric radii. 
HNC behaves oppositely, peaking in the 40 km s 1 feature which 
is the strongest individual cloud feature seen along the line of 
sight. In this regard, CF + most closely resembles HNC. 

In Fig. 3 we compared the column densities of CF + , HCO + , 
C 2 H and HNC more quantitatively along all the sightlines where 
CF + has been observed in absorption. The 33 kms -1 cloud to¬ 
ward W49N has molecular column densities most nearly similar 
to the sightline toward BL Lac at Ay = 1 mag, but only 40% as 
much CF + . The 39 km s~ 3 feature toward W49N is similar to the 
higher-column density translucent sightline toward 3C 111 (Ay 
<; 4 mag), including CF + . N(CF + ) shows a quadratic dependence 
on N(C 2 H) and most closely resembles HNC in its detailed be¬ 
haviour. 

As shown in Table 4 the relative abundance X(CF + ) = 
N(CF + )/N(H 2 ) observed in absorption toward W49N are 3.3 x 
10' 10 for the gas around 39 kms' 1 and 1.4 + 0.16 x 10' 10 for 
the other four velocity intervals in which CF + was detected. To 
interpret these results we embedded the fluorine chemistry in a 
model of the heating, cooling and H 2 and CO formation in dif¬ 
fuse and translucent gas using a recent remeasurement of the F 
+ H 2 —> HF + H formation reaction. Using CH as a surrogate 
for H 2 with X(CH) = 3.5 x 10 3 the models correctly reproduce 
the observed column density and relative abundance of HF with 
X(HF) = N(HF)/N(H 2 ) = 1.1 x 10 8 . However, they also predict 
N(CF + )/N(H 2 ) = 5.3 x 10~ 10 and N(CF + )/N(HF) ~ 1/20 which 
is too large by factors of four for CF + seen over most of the ve¬ 
locity range seen toward W49, or by 40% for the feature around 
39 kms- 1 . 

Given the correctness of the predictions for HF, we hypoth¬ 
esized that the discrepancy in the predicted CF + abundance and 
CF + /HF ratio in the models results from the absence of a sig¬ 
nificant CF + destruction process. This process is presumably 
suppressed when the CF + relative abundance is unusually high. 
We surmised that photodissociation of CF + is important and we 
showed, by examination of the recombination rates in the mod¬ 
els, that a photodissociation rate F ^ 10 9 s -1 (comparable to 
that of CH) could cause a substantial decrease in the CF + abun¬ 
dance, introducing number and column density dependences 
that otherwise are largely absent. However, the free space pho- 
todissocation rate has recently been calculated and found to be 
T = 2 x 1 O ' 10 s _1 , far too small to explain the observed dis- 
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crepancies. In any case, uncertainties remain in the adopted CF + 
formation and recombination rates, which may need revision. 

In Fig. 4 we compared the optical depth profiles of C + and 
HF, the species that form CF + , with that of CF + itself. The pro¬ 
file of r(C + ) varies little across the peaks in r(HF) and r(CF + ), 
owing to blending of atomic and molecular gas and somewhat 
poorer velocity resolution, with, perhaps, a decreased C + frac¬ 
tion in the high-column density cloud. If C + is the dominant form 
of gas-phase carbon it dominates both the formation and de¬ 
struction of CF + , leaving CF + to trace N(HF) and NtFL). In this 
experiment along the complex and heavily-blended low-latitude 
line of sight to W49N, CF + signaled the coexistence of C + and 
HF, but primarily traced the peaks of the PL distribution. 


Appendix A: Extraction of the CF + line profile 

As noted in Sect. 2.1, We extracted the CF + profile in the pres¬ 
ence of an overlying H66e recombination line. We have confi¬ 
dence in the identification of this extraneous feature by compar¬ 
ison with maps we made of the more isolated H53y line, whose 
38 kms" 1 linewidth is the same as that of the gaussian profile 
fitted to extract the CF + profile. 

Shown in Fig A. 1 is the observed continuum-subtracted pro¬ 
file, as derived from a linear fit to the spectral baseline well out¬ 
side the velocity interval shown. A gaussian profile with FWHM 
38.2+1.6 km s~' was fit to the spectrum using spectral channels 
within velocity intervals chosen to coincide with regions of weak 
HF absorption, in order to avoid removing any broad underlying 
component of CF + absorption. The fitted gaussian is shown in 
red and the velocity intervals used for fitting are outlined as blue 
dashed rectangles. The sum of the gaussian and linear fits was 
then used to derive the final continuum-subtracted profile from 
the original data. The continuum flux at the observed position 
was 3.47 Jy and the rms line/continuum noise at zero-absorption 
in the final continuum-subtracted spectrum was 2.4 10 5 . 

We tested the reliability of the fit by adding to the data vary¬ 
ing amounts of absorptio n mimicking the optical depth spectrum 
of HF (iGerin et all 120I2L to which CF + is closely related, and 
then observing whether the added absorption was present in the 
derived CF + profile after gaussian decomposition. This was the 
case as long as the velocity intervals used in the profile decom¬ 
position were those in which HF absorption toward W49N is 
weak. We are confident that no broad component of CF + absorp¬ 
tion was removed by the profile decomposition, although such a 
feature may be well be present below the levels of detectability 
of the experiment. 


Appendix B: Line of sight kinematics 

The line of sight kinematics toward W49 are illustrated in very 
compact form in Fig. B.l for a flat rotation curve, as noted in 
the figure caption. W49N is at the far kinematic distance only 
slightly inside the Solar Circle so the same path is sampled in 
both emission and absorption. The maximum velocities of fea¬ 
tures observed in emission and absorption correspond well to 
that expected from the model velocity field. The so-called ter¬ 
minal velocity feature at v ^ 60 l< m s “ 1 occurs as material over 
a wide range of line of sight distance collects in a narrowed ve¬ 
locity range owing to the smaller line of sight velocity gradient. 
The velocity range up to +20 kms -1 includes the nearest 1.5 
kpc along the line of sight and constrains the absorption to arise 
inward of 7.5 kpc in galactocentric distance.. 



Fig.A.l. CF + and recombination line profiles observed toward 
W49N. The solid black histogram is the continuum-subtracted 
profile using a linear spectral baseline fit to signal-free channels 
well outside the region shown. The red line shows the gaussian 
profile fitted to the recombination line emission over the spec¬ 
trum segments outlined in blue, which were chosen to coincide 
with regions of small HF, hence CF + absorption. 
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Table 2. Optical depth integrals over selected velocity intervals 11 


Vel. range 
km s 1 

c 3 h 2 

km s -1 

HCO + 
km s -1 

HNC 

km 

red CF+ 
km s 1 

HCO 
km s 1 

H l1 CO + 

km 

CH 
km s -1 

23.9..31.3 

0.520(0.040) 

— 

0.000(0.017) 

0.0000(0.0027) 

0.0000(0.0121) 

— 

0.035(0.019) 

31.6..35.7 

1.380(0.040) 

1.950(0.018) 

0.270(0.018) 

0.0110(0.0024) 

0.0540(0.0096) 

0.0000(0.0112) 

1.144(0.014) 

35.7..42.0 

3.140(0.070) 

6.210(0.074) 

1.970(0.036) 

0.0514(0.0029) 

0.1322(0.0115) 

0.1240(0.0140) 

2.040(0.019) 

42.1..50.0 

0.660(0.040) 

0.840(0.016) 

0.000(0.023) 

0.0000(0.0032) 

0.0325(0.0126) 

0.0000(0.0154) 

0.861(0.017) 

50.0..56.8 

1.570(0.045) 

3.980(0.025) 

0.400(0.022) 

0.0161(0.0030) 

0.0867(0.0118) 

0.0000(0.0144) 

1.660(0.017) 

56.8..61.5 

2.780(0.054) 

7.270(0.070) 

1.200(0.023) 

0.0255(0.0026) 

0.2410(0.0102) 

0.1234(0.0123) 

2.150(0.018) 

61.5..66.1 

2.370(0.053) 

4.490(0.036) 

0.880(0.021) 

0.0218(0.0025) 

0.1410(0.0097) 

0.0649(0.0118) 

1.830(0.017) 

66.1..70.3 

0.628(0.330) 

0.810(0.013) 

0.110(0.017) 

0.0000(0.0024) 

0.0275(0.0095) 

0.0000(0.0113) 

0.670(0.013) 


“ Quantities in parenthesis are lcr statistical error estimates, see Sect. 2 


Table 3. Line frequencies and opacity-column density conversion factors 



c 3 h 2 

HCO + 

HNC 

CF + 

HCO 

H 18 CO + 

C 2 H 

HF 

C + 

CH 

frequency" 

factor* 

85.339 

4.36E12 

89.189 

1.12E12 

90.664 

1.76E12 

102.588 

1.35E13 

86.671 

2.26E13 

86.754 

1.13E12 

87.317 

6.52E13 

1232.476 

2.404E12 

1900.537 

1.40E17 

532.774 

3.64E13 


“ GHz 

b units for all entries are cm~ 2 (km s -1 ) -1 


Table 4. Relative Abundances" 


Vel. range 
km s' 1 

C 3 H 2 

HCO+ 

HNC 

CF+ 

HCO 

C 2 H 

HF 

31.6..35.7 

1.48E-09 

1.84E-09 

3.99E-10 

1.25E-10 

5.08E-11 

3.99E-08 

1.30E-08 

35.7..42.1 

1.66E-09 

3.28E-09 

1.63E-09 

3.26E-10 

6.98E-11 

4.56E-08 

1.29E-08 

42.0..50.0 

1.03E-09 

1.05E-09 

<1.36E-10 

<1.44E-10 

4.07E-11 

2.89E-08 

6.77E-09 

50.0..56.8 

1.23E-09 

2.58E-09 

4.08E-10 

1.26E-10 

5.62E-11 

2.97E-08 

9.57E-09 

56.8..61.5 

1.45E-09 

3.64E-09 

9.45E-10 

1.54E-10 

1.21E-10 

3.60E-08 

1.20E-08 

61.5..66.1 

1.43E-09 

2.64E-09 

8.14E-10 

1.54E-10 

8.30E-11 

3.17E-08 

9.68E-09 

66.1..70.3 

1.25E-09 

1.30E-09 

2.78E-10 

<1.39E-10 

4.42E-11 

2.66E-08 

8.59E-09 

Suit? 

1.41E-09 

2.66E-09 

7.89E-10 

1.57E-10 

7.43E-11 

3.54E-08 

1.08E-08 

Mean" 

1.36E-09(0.21) 

2.33E-09(0.98) 

7.46E-10(5.05) 

1.77E-10(0.85) 

6.65E-11(2.82) 

3.41E-08(0.68) 

1.04E-08(0.24) 


“ Quantities tabulated are N(X)/N(H 2 ) where N(H 2 ) = N(CH)/3.5 1CT 8 
b Entries labelled Sum are ratios of the total column density over all intervals to the total CH 
c Entries labelled Mean are the average(rms) of the detections in the body of the table. 
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Fig.B.l. W49N line of sight kinematics in a nutshell. The mot¬ 
tled blue line shows the galactocentric radius R at each velocity 
assuming a flat rotation curve with 0(7?) = 0(7?o) = 220km s ', 
R(> = 8.5 kpc. The solid green curve shows the line of sight dis¬ 
tance or distance s p at ea ch velocity and the distance to W49 p, 
11.1+0.8 kpc f Zhang et al.L 1 2013 1 is illustrated by an an asterisk. 
The dotted red curve shows the line of sight velocity gradient 
dV/dp due to Galactic rotation. Shown superposed to illustrate 
the observed range of spectral features are the CO profile in units 
of Kelvins, offset but not scaled, and the CF + absorption spec¬ 
trum scaled upward a factor of 100 and offset, both from Fig. 
1 . 
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